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Abstract 

/3-functions for abelian and non-abelian gauge theories are studied in the regime 
where the large N flavor expansion is applicable. The first nontrivial order in the 1/N 
expansion is known for any value of Na, and there are also various indications as to 
the nature of higher order effects. The singularity structure as a function of Na has 
implications for the existence of nontrivial fixed points. 

For a sufficiently large number of flavors a non-abelian gauge theory will loose asymptotic 
freedom and will in this way resemble an abelian gauge theory. We shall thus focus our 
attention on the abelian case and then later extend the discussion to the nonabelian case. 

It is generally believed that a U(l) gauge theory with iV charged fermions has a running 
coupling that grows monotonically towards the ultraviolet, and thus suffers from a Landau 
pole. This is the indication from the one-loop /3-function. But there is much more known 
about the perturbative /3-function and there have been recent calculations that have extended 
our knowledge to 5-loops. We also have complete knowledge of the first nontrivial order in 
the 1/N expansion for any value of Na, and so this makes the large N expansion a useful 
way to organize the perturbative expansion. The question is whether any of this allows us 
to glean anything further about the possible existence of nontrivial fixed points. Although 
any perturbative approach will introduce a renormalization scheme dependence, it can still 
be hoped that the existence or nonexistence of a fixed point will leave some mark on the 
perturbative results. 

According to a lattice result pQ there is no nontrivial fixed point in a £7(1) gauge theory 
for N = 4. We shall be concerned with larger N where the large N expansion suggests other 
possibilities. This provides some motivation to study larger N values on the lattice as well, 
for sufficiently large values of Na. Extending the current lattice result to N = 8, 12, 16,... 
would appear to be relatively straightforward. 
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The one loop result is /3(a) = 2A/3 where A = Na/n. We may write an expansion in 1/N as 
follows, 

3 /3(a) , , ^F,(A) 
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The "1" corresponds to the one loop result and we shall refer to it as the zeroth order term in 
the 1/N expansion. Each Fi(A) represents a class of diagrams having the same dependence 
on iV when A is held fixed, and such diagrams exist to all orders in A. If the functions 
were bounded then for sufficiently large N one could conclude that the zeroth order term 
dominates and that the Landau pole is unavoidable. But singularities in the Fi(A) will keep 
us from drawing this conclusion. 

We collect together what is known about the Fi(A)'s in the MS renormalization scheme. 

Fi(A) = / h(x)dx (3) 
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Important for our study is the fact that Fi(A) is known completely [2j. We have expressed the 
integrand I±(x) in a form that makes more clear the location of its zeros and poles. The A 3 
terms in F 2 (A) and F^(A) were calculated in [3], the A A term in F 2 (A) in [4], and the ^(^4) 
term in [5]. The latter two results are 5- loop calculations. 

One way to express the results in ([3]{7]) is to plot their sum and ignore what is not known. 
The result for the 3(3(a)/2A for various N is displayed in Fig. (1). A zero would indicate a 
nontrivial fixed point, but the zeros are occurring at values of A that are too high to ignore 
higher order terms. Thus we cannot deduce much from this plot, except to notice sensitivity 
of the /3-function to N. 

The 2- loop contribution to f3(a) involves one fermion loop and one internal photon and 
it gives rise to the first term in the expansion of Fi(A), which is |A The higher order 
terms in F\ (A) correspond to the insertion of the appropriate number of fermion loops on the 
photon line, and these bubble chains have been summed up to produce the result ^ [2J. An 
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important feature of Fi(A) is that its expansion displays a finite radius of convergence, which 
is nonvanishing due to the slower than factorial growth in the number of diagrams. We expect 
this to be true of the other Fi(A)'s as well. 

Simple poles of alternating sign appear in I\{x) at x — | + n for integer n > 0. The 
integration can be handled with a Cauchy principal value prescription and the result is shown 
in Fig. (2). Clearly Fi(A) only changes logarithmically as the singular points are approached. 
And these singularities become weaker for larger n. Close to the first singularity at A = 15/2 
we find 



Fi(A) « — »(log(l - —A)) + 0.3056. 



Nevertheless even this weak singularity can cause the /3-function at this order to vanish at a 
fixed value of N. As Fig. (2) indicates there will be two nearly coincident zeros of l + Fi(A)/N 
at 

A= y iO.OUTe- 15 ^ 7 . (9) 

The lower (upper) one is an ultraviolet (infrared) fixed point. In either case the running 
coupling achieves its fixed point value at some finite scale /z*, and at this scale the running of 
the coupling abruptly stops. 

It is useful to compare the /3-function to the 7 m -function defined as 



t»w = -tot < 10 ) 

Gi(A) is also known to all orders in A, and in fact it is directly related to Fi(A). From the 
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Figure 2: -Fi(A) as defined in 



results of [2] one can deduce that 

dFx(A) 
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Thus the singularities in these two functions occur at the same locations. This could be 
expected since the bubble chain re-summation is intrinsic to both functions. The strength of 
the singularities are also related; the logarithmic singularities in F\(A) correspond to simple 
poles in G\(A). 

It could be argued that the Cauchy principal value prescription used in the evaluation of 
Fi(A) is not unique. Rather than approaching the pole equally closely from the two sides, 
another prescription would be to approach the pole unequally from the two sides. This would 
shift the /3-function on the right of a singularity by an additive constant, as allowed by ( JlT| . 
But this ambiguity does not alter the appearance of the fixed points at first order in 1/JV. 

It is important to know how large N has to be for the 1/N expansion to be under control. 
We can require that the known expansion terms of the higher Fi(A)'s be sufficiently small for 
A as large as the radius of convergence of Fi(A), at A = 15/2. If we rescale A = 15/2A and 
N = 16 N then the expansion &2\) numerically reads 



1 + ^(.3516i-.8057i 2 - 1.567l 3 + 5.342i 4 + 1.60i 5 -15.9i 6 + ...) 
N 

+ i(-.0206i 2 - 1.602i 3 - 3.244A 4 + ...) 
N 2 

+ i(-.0555i 3 + ...) + i(.1198i 4 + ...) 
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We see that N > 1 or more is needed for the terms of successively higher powers of 1 / N to be 
under control. One can in particular compare the leading terms at each order in 1/N, i.e. the 
A 1 term at order 1/N 1 . These terms correspond to the one fermion loop diagrams, and they 
have a special significance in that they are renormalization scheme independent [3]. 

Thus for sufficiently large N the higher orders in 1/N are under control in the usual 
sense of an asymptotic series. But the presence of singularities in the Fi(A) indicates that 
the 1/N expansion needs to be reconsidered for A close to these singularities. Although the 
singularities of F\{A) are logarithmic, we do not expect this to continue for the higher Fi(A). 
The appearance of poles in G\(A) reinforces this view. In particular if F%(A) has a simple 
pole at A = 15/2 then this should completely dominate the logarithmic singularity of F\(A), 
as the exponentially small spacing in ^ makes clear. 

It is interesting that the first three coefficients in the expansion of F 2 (A) are negative. 
This is certainly consistent with a pole, since the expansion of the expression 

^OTTli + OMO 
A - 0.494 

has the same first three terms. But we would argue that this does little to prove the existence 
of a pole. In the Fz(A) expansion we note that the 2 fermion loop A 3 term in F 2 (A) receives 
contributions from graphs with topology different from the 1 fermion loop A 2 term. That is, 
unlike the case of Fi(A), dressing photon lines of the A 2 graphs does not give all A 3 graphs. 
In fact the graphs with the new topology (the light-by-light scattering contributions) give 
the dominant contribution to the A 3 term [3]. In this situation we do not expect a Pade 
approximant to be very predictive until more orders in the expansion are known. 

An example of the appearance of a pole in an exact /3-function is provided by SU(N C ) 
SUSY pure gluodynamics where P Ej 

n, x 3a 2 1 a c .„ 

m = ~iri=w? (14) 



A pole of this sign, as in (13), means that the coupling will evolve towards the singularity in 
the infrared, and it does so from both the weak and strong coupling sides of the pole. The 
coupling reaches the singularity at some minimum finite renormalization scale and the 
gauge theory ceases to provide a description below this scale. This suggests that the theory 
develops a mass gap and/or some fundamentally different description is needed for energies 
below [A*. The authors of [8 J argue that there is evidence of such a pole in the QCD (3- 
function from the study of Pade approximants. A pole of opposite sign would instead produce 
an ultraviolet cutoff on the gauge theory description. 

Beyond the A 3 term in ^(A), including the known A A term, the diagrams are obtained 
by simply inserting fermion loops into photon lines, thus building up bubble chains. This is 
true for sufficiently high powers of A in the expansion of any Fi(A); one ends up only dressing 
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photons lines with fermion loops in graphs that belong to a basic set of topologies. It is for 
this reason that all these functions should have a nonvanishing radius of convergence, and the 
bubble chains should generate singularities in the higher Fi(A) just as they did for F\(A) and 
G\(A). We expect that the set of locations of singularities in the higher Fi(A) will include 
the locations of singularities in Fi(A). 

Let us consider a set of poles in F 2 (A) occurring at the same locations at the singularities 
in F\(A) (we shall provide more evidence for this shortly). Poles in F 2 (A) at A n = y + 3n 
need to be considered along with the zeroth order term in (pi). Depending on the sign of these 
poles, the result is the existence of various nontrivial infrared and/or ultraviolet fixed points. 
For example an infrared pole at A = 15/2, i.e. having the sign of 1/(A — A Q ) as in (13) and 
(14), implies an ultraviolet fixed point at some Auv < A (see Fig. (3a)). When the coupling 
is below Auv the theory flows to the known weakly coupled behavior in the infrared, but in 
the ultraviolet the Landau pole has been eliminated in favor of a fixed point. Possibilities for 
nontrivial fixed points in the regions between successive poles at A^ and A i+ i are shown in 
Figs. (3b,c,d). 

We now describe some evidence that is relevant not only to the existence of poles in F 2 (A), 
but also to their location and sign. The first input comes from the 7 m -function, since it turns 
out that G 2 (A) is also known to all orders in A. This comes from an impressive calculation 
in QCD |9j, from which the corresponding QED result can be extracted. The surprisingly 
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simple result is that the singularities in G 2 (A) occur at the same locations as in Gi(A), and 
each simple pole in G\{A) has been replaced by a double pole in G 2 (A). 

Further input comes from examples of theories where the analogs of G 2 (A) and F^{A) 
have the same singularity structure; that is the set of poles of these functions have the same 
location, sign and order of pole. In particular in N flavor (<ft 2 ) 2 theory in d = 4 — 2e dimensions 
the critical exponents X(e) and ui(e) have a correspondence to 1 — 7 m and 2(3'. In jlO] these 
critical exponents are encoded in the corresponding critical exponents of the large N a-model 
in d = 2 — 2e dimensions, where the functions appearing at order 1/N 1 are labeled Aj(e) and 
u)i(e). These functions are the analogs of — Gi(A) and F[(A). From the results in [10] the 
functions — X 2 (e) and lu 2 (e) have singularities that match in location, sign and order of pole, 
as do —Xi(e) and Wx(£:)q 

Thus it would not be surprising for QED to display similar behavior, so that G2{A) and 
F^(A) exhibit a common singularity structure just as G\(A) and F[(A) do. Such behavior 
could be expected due to the simple relationship between the graphs contributing to Fi(A) and 
Gi(A), suggesting that the behavior also extends to higher orders. Cutting out an external 
gauge field vertex from any graph contributing to Fi(A) (vacuum polarization graphs) and 
replacing the other external vertex with a mass insertion gives a graph contributing to Gi(A) 
(mass renormalization graphs). This is not true in QCD due to the gluon self- interactions, 
and so there is no corresponding behavior for QCD even at first order in 1/N. 

Thus supposing that F^{A) and G 2 (A) have the same singularity structure, the signs of 
the poles in F 2 (A) at A n = y + 3n are determined by the known results [9] for G 2 (A). This 
would imply that the poles in F 2 (A) all have the same sign and that they are all ultraviolet 
poles, i.e. opposite in sign to (13) and (14). Thus of the various possibilities displayed in 
Fig. (3), the available evidence suggests that only Fig. (3d) is correct for the interval between 
any two adjacent singularities. 

If this is correct then at this order in the 1/N expansion there is an infinite number of 
new theories each flowing towards an associated nontrivial infrared fixed point. There is 
still the weakly coupled theory that flows to the trivial infrared fixed point. All the theories 
are cut off in the ultraviolet by the UV poles. We also note that for large N some of the 
nontrivial infrared fixed points can occur at values of A smaller than the critical value for 
chiral symmetry breaking, A CT \ t = N/3. 

To proceed to higher order we must once again reconsider the 1/N expansion close to the 
singularities. As a singularity is approached the /3-function experiences a significant change, 
and a change of 3(3/2 A of order unity would be an indication that the 1/N expansion is 
breaking down for these values of A. The higher Fi(A)'s must be considered in these regions. 
For these functions to overcome the suppressions from large N and loop factors, they would 



1 wi(e) has simple poles, while w 2 (e) has an additional pole and all its poles are fourth order. 
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need to have progressively higher order poles. 

In fact it is plausible that the Fi(A) function has (i — l)th order poles, since this behavior 
is already suggested by the increasing order of the poles in the G\{A) and G 2 {A) functions. 
A 1/N expansion with poles of ever increasing order would indicate the existence of a set of 
essential singularities in the /3-function. But the implications for fixed points remain similar to 
our discussion of simple poles. Whenever an essential singularity drives 3/3/2A by an amount 
of order unity or more in the right direction then a nontrivial fixed point can result. 

We may speculate further about the nature of the higher order poles based on an emerging 
pattern for poles of even or odd order. The adjacent logarithmic (even order) singularities 
in Fi(A) are alternating in sign while the simple (odd order) poles in F 2 (A) are all of the 
same sign. Let us consider the continuation of this pattern (alternating even order poles and 
same sign odd order poles) to higher orders, where we treat the overall sign at any order as 
unknown. It means that on a given interval between adjacent singularities, and at any given 
order, 3/5/2A diverges with opposite signs on the two ends of the interval. The divergences 
of the all-orders-summed result (assuming that the summed result diverges) can then also 
be expected to come with opposite signs on the two ends. This would imply at least one 
nontrivial ultraviolet or infrared fixed point on each such interval. Also, the sum of the odd 
order poles produces a divergence pattern that is the same on every interval, while for the sum 
of the even order poles the pattern flips in sign on adjacent intervals. Thus the existence of 
infrared fixed points in the total result, at least on every second interval, is made more likely 
if the odd poles are always of the ultraviolet variety as in Fig. (3d). Due to cancellations it 
is also conceivable that 3f3/2A approaches finite or even vanishing values at the endpoints of 
some intervals. Such behavior that depends on the side of approach can also be consistent 
with essential singularities. 

We now turn to a SU(N C ) gauge theory, where results are expressed in terms of the 
Casimirs, C G = N c , C R = (iV c 2 - l)/2N c , T R =1/2. N still denotes the number of flavors. 
The term in the /3-function of zeroth order in the 1/N expansion is as before, (3(a) = 2A/3 + ..., 
if we choose a new definition for A = NT R a/ir. Similarly the 1/N expansion is 

3 /3(a) J^Hi(A) 

2^- = 1 + L^- ( 15 ) 

i=l 

From the results in [TT] we can deduce tha 



11 c r A/3 

H\(A) = ---^+/ h{x)h{x)dx, (16) 



4 T R Jo 

j / \ C R (20-43x + 32x 2 -Ux 3 + 4x 4 ) C g 
2[X) T R + 4(2*-l)(2*-3)(l-*2) T R [U) 



2 This result is presented somewhat more explicitly than in [llj . but its series expansion is in agreement. 
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Figure 4: Hi(A) as defined in (16) for N c = 3. 



h(x) is defined in Q and so up to the definition of A and the Cr/T r factor, the Cr term is 
just the QED result for Fi(A). By inspection of h(x) and h(x) one can see that the Cq term 
brings in a new pole in the integrand at x = 1 (A = 3). This pole is of the same sign as the 
pole at A = 15/2. We show a plot of Hi(A) in Fig. (4) for N c = 3. Compared to the Fx(A) 
function of QED, Hi(A) is negative and its first singularity is occurring at a smaller value of 
A (and a). 

As in the QED case, a key question is how large N has to be for the 1/N expansion to be 
under control, for A at least as large as 3. It is apparent that higher N will be needed, since 
factors of iV must now compete with factors of N c . If we rescale A = 3A and N = 32N then 
the expansion (15) for N c = 3 numerically reads 



1 + 



+ 



+ 



N 
1 

AT3 V 



-0.5156 + 0.8906i - .6348i 2 - .2372A 3 + .4372i 4 - 0.0994i 5 - .0912i 6 + ...) 



'-0.224L4 + 0.9216i 2 - 2.003A 3 + 



.1471A 2 + 1.073A 3 + 



1 , 



.1412A 3 + 



The leading terms of the expansions of H 2 (A), H 3 (A), H±(A) have been obtained from the 
4-loop results in [12J. Due to the pure glue contributions, each order begins with a lower power 
of A as compared to QED. The main observation here is that large values of N are needed, 
iV > 32 or more, for the 1/N expansion to be under control. 

For such large values of N, asymptotic freedom has been lost, and the implications of 
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Figure 5: 3(3(a)/2A at first order in l/N where a zero occurs at the critical coupling for chiral 
symmetry breaking (N = 8.88 and N c = 3). 



the singularities in H\(A) and the higher Hi(A)'s will resemble that of QED. The large N 
required also indicates that the large N expansion lacks quantitative control for the study 
of other interesting phenomena that occur for values of N 32]^] One such question is 
the lowest value of N at which the Banks-Zaks fixed point survives before chiral symmetry 
breaking occurs. 

For completeness we give the answer to this question while only keeping the 1 + H\(A)/N 
terms in (15). We can find the N such that the zero of 1 + H\(A)/N occurs at the usual 
estimate for the critical coupling for chiral symmetry breaking, or A crit = NTr/3Cr. For 
N c = 3 this gives N = 8.88 as the lower bound on the conformal window. This may be 
compared with the conventional result of iV = 11.9 which is derived from the 1- plus 2- loop 
terms of /3(a). We plot 1 + H^A)/N for N = 8.88 in Fig. (5). 

At order l/N 2 , the singularity structure of the QCD /3-function is significantly more com- 
plicated. And unlike QED, the singularity structures of the 7 m and /3-functions do not match. 
For example H[(A) has the additional singularity at A = 3 when compared to the QCD 7 m - 
function at order l/N, which is essentially the QED result Gi(A). As mentioned before, this 
can be understood diagrammatically in terms of additional contributions to the /3-function 
due to the gluon self-coupling. Nevertheless there will be contributions to singularities in 
H' 2 {A) that can be associated with the singularities that appear at order l/N 2 in the QCD 
7 m -function. Since the latter is known [9] we can deduce that the H2(A) singularity struc- 



3 Of course if the study involved values of A significantly less than ^4 = 3 then such large values of N would 
not be needed. 
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ture will contain the following contributions: a simple UV pole at A = 3/2, a positive log 
singularity at A = 3j^] a negative log singularity at A = 9/2, and IR poles of third order at 
A = 15/2 + n. Some of these singularities could change sign or be replaced by even stronger 
singularities in the full result for H2(A). The UV pole is interesting since it would imply an 
infrared fixed point just above A = 3/2. 

In summary we have considered the implications of singularities appearing in the large 
N flavor expansion of the /3-functions of QED and QCD. The logarithmic singularities that 
appear at first order in 1/N are the first signs of simple and higher order poles that will 
appear at higher orders in 1/N. In the case of QED these singularities are expected to be 
related to those of the 7 m -function, for which the singularity structure is surprisingly simple 
at order 1/N 2 . This gives information about the signs and locations of a set of poles in the 
/3-function at this order. This in turn implies an infinite set of nontrivial infrared fixed points. 
The poles are likely to turn into essential singularities when all orders of the 1/N expansion 
are considered, but the existence of nontrivial fixed points will persist. For QCD we have 
indicated an even richer singularity structure for which we have less information. 
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